A Predictive Minimal Model for Neutrino Masses and Mixings 
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A model is considered in which the scale of the heavy singlet neutrinos is a few orders of magnitude 
below the grand unification scale and where right-handed vector bosons play still a negligible role. 
In a basis with diagonal up-quark and Dirac-neutrino mass matrices it is assumed that the heavy 
neutrino mass matrix has only zero elements in its diagonal, in analogy to the light neutrino mass 
matrix in the Zee model. Connecting then the remaining matrix elements with the small parameter 
describing the hierarchy of quark masses and mixings and by assuming commutativity of the charged 
lepton with the down-quark mass matrix, the calculation of all neutrino properties can be performed 
in terms of the two mass differences relevant for atmospheric and solar neutrino oscillations. CP- 
violation is directly related to CP-violation in the quark sector. 



I. INTRODUCTION 



Recent experiments give strong evidence that neutrinos oscillate and have finite masses [§J , |§] ■ These experiments 
provide the first conclusive evidence for physics beyond the standard model. Good candidates for a corresponding 
extension of the standard model are supersymmetric or non-supersymmetric grand unified theories For recent 
summaries in which many models are discussed and for further literature see Q. Detailed grand unification models 
need many assumptions about numerous new particles and their representations, Higgs potentials and couplings. One 
, obtains interesting relations, but it is often hard to see to what extent the many assumptions are really necessary 
and decisive for the light neutrino sector. In this article I will formulate a model in which the heavy neutrinos get 
\Q , masses several orders of magnitude below the scale of grand unification. The possibility exists that these masses are 
C*"** ■ generated dynamically in some analogy to the generation of light neutrino masses in the Zee model J5|. Only few 
particle properties at the unification scale are then directly relevant for the neutrino problem. The minimal model 
obtained in this way provides an example in which the significance of the assumptions are obvious and the neutrino 
properties can be calculated in all details. A disadvantage compared to an explicit treatment of a specific grand 
unified theory is that there are no immediate additional predictions about new physics at the unification scale. 
Oscillations involving 3 light neutrinos^] can be described in terms of 6 parameters: 2 parameters for the difference 
'■^ \ between the square of the neutrino masses and 4 parameters for the unitary matrix U which relates - in a basis in 
which the charged lepton mass matrix is diagonal - the flavor eigenstates v e ,i/^, v T with the mass eigenstates v\,i>2, v% 

q.: 

<D . v a = s £ j XJ al v t (1) 



The theoretical understanding of neutrino physics and of these parameters in particular requires, however, the con- 
sideration of many more - a priori unknown - quantities. An appealing way to extent the standard model is to add 
to the 3 two-component fields new two-component fields v e , , v r , which are singlets with respect to the standard 
model gauge group and acquire large masses corresponding to the scale of the new physics involved. The well-known 
see-saw mechanism then provides for the observed smallness of the light neutrino masses Q . 

There are then 3 mass matrices to consider]^]: the Dirac neutrino mass matrix m® lrac which connects the old with the 
new two-component fields, the mass matrix for the charged leptons ttle and the mass matrix of the singlet neutrinos 

Mr. 

There is some freedom for choosing a basis for the 3 matrices: We may take the Dirac neutrino mass matrix diagonal 
and real which leaves us with 3 parameters for this matrix. Below the scale of new physics the charged lepton mass 
matrix can be taken hermitian. After a proper adjustment of the phases of the lepton fields this matrix contains 7 
parameters. Finally, the symmetric mass matrix Mr can be transformed to have real diagonal elements and contains, 
therefore, 9 parameters. 



* E-mail: B.Stech@thphys.uni-heidelberg.de 

x In this article the possibility of the existence of light sterile neutrinos is not discussed even though they are needed in a 
general analysis, in particular, if the finding of the LNSD group || are included. 

2 A classification of three-neutrino models dealing with these mass matrices is provided in ref. M. Ansatze related to the 
present approach can be found in ref. |8j. I will closely follow and extend my previous publications [M on this subject. 
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In the next sections we will restrict these 19 parameters using the ideas mentioned above and by suggesting 
connections with the known structure of the mass matrices of up and down quarks. The see-saw formula 

my = -m° trac ■ M^ 1 ■ (m^ rac ) T (2) 

together with the mass matrix for charged leptons then allows to calculate all the properties of light neutrinos of 
interest here. 



II. THE DIRAC NEUTRINO AND THE UP-QUARK MASS MATRICES 

Below the scale of new physics we are free to choose a basis in which the Dirac neutrino mass matrix m® lrac and the 
up-quark mass matrix mjj are diagonal simultaneously. In grand unified theories (taking low Higgs representations) 
one finds that at the unification scale m\j and m® lrac are closely related and could even be equal there jf(|. We 
will assume, therefore, that at the mass scale of the heavy neutrinos, which will turn out to be much lower than the 
unification scale, m^ lrac and mjj can be diagonalized simultaneously. In section 4 we will show that this mass scale 
which we denote by Mo is of order 10 11 GeV. 

The observed hierarchical structure of charged fermion masses and mixings allows to express the corresponding mass 
matrices in terms of powers of a small quantity. Using for this parameter a = (nic/nit) 1 ' 2 , the diagonal up-quark 
mass matrix can be written |j] 

fa 4 \ 
my (M ) = o 2 m t (M ) 
V ij 

a = 0.057 . (3) 

With the help of the renormalization group equations Eq. (|^) together with the down-quark and charged lepton mass 
matrices (see section 5) provides values for the masses of up, charm and top quarks at the scale of the vector boson 
Z which agree with the known result fIH| 

m u {m z ) = 1.9 ± 0.4 MeV, m c {m z ) = 0.61 ± 0.05 GeV, 

m t {m z ) = 173 ±5 GeV, 
m d {m z ) = 3.4 ± 0.6 MeV, m s {m z ) = 0.064 ± 12 GeV, 

m b (m z ) = 2.90 ± 0.04 GeV (4) 

within the given error limits. For m® lrac at the scale Mq we can write therefore 

ycr 4 

m^ rac (M ) =| xa 2 \ m t (M ) (5) 
1 

with x = 0(1), y = 0(1). There is no point in taking an additional O(l) parameter for the 33 element. As can be 
seen below this would only affect the scale parameter Mq which is a fit parameter in our model. 



III. THE MASS MATRIX FOR THE HEAVY SINGLET 
NEUTRINOS 



In the absence of a detailed theory describing new physics beyond the standard model, one can only speculate about 
the mass matrix for the heavy neutrinos. However, one can hope that the small parameter a which governs the up and 
down quark mass matrices plays also here an important role j^] . Because of the self-coupling of the heavy neutrinos 
the assignment of generation quantum numbers to those fields ||l2|| gives additional restrictions for the powers of a 
occurring in this matrix. As an Ansatz for Mr I will use here a matrix with only zero elements in its diagonal (see 
also fl3|| ), in analogy to the matrix for the light neutrinos obtained in the Zee model g. Mr could have a similar 
origin, this time, however, involving the "right-handed" fields in a right-left symmetric theory such as 5O(10) or E6 
at the grand unification scale S> Mq- The elements of a matrix of this form can easily be connected with the elements 
of the mass matrix m® %rac (and mu) with the powers of a related to appropriate generation charges of the singlet 
neutrino fields. Therefore, I propose 
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Mr = xa A a M (6) 



a 



with a = 0(1). In this Ansatz, the ratio between (Mr)\% and (Mr) 13, namely xa 2 , has been taken to be identical 
to the ratio between the elements (m® trac )i2 and (m^ trac )z7i in order to have a close relation with m® lrac . Since 
this ratio is a real number, all elements of Mr can be taken to be real by a proper phase choice for the singlet 
neutrino fields. The powers of a occurring in (^|) correspond to generation charges for the fields i> e , v T equal to 
5/2,3/2,-1/2, respectively, when setting the generation charge of the scalar field which produces the heavy masses 
equal to —1. 

Mr can be diagonalized by an orthogonal matrix Vr defined by 

Mr = V R M d R ia9 ° nal Vl . (7) 

The two large eigenvalues of Mr differ up to order a 2 only in sign (like the eigenvalues of a Dirac neutrino mass 
matrix) (~ ±aMo) while the third eigenvalue is smaller by the factor er 4 . Because of the almost degeneracy of two 
eigenvalues it is interesting to look at the up-quark mass matrix when taken in a basis in which Mr is diagonal. One 
finds - up to order a 2 - for x = 1 



VrttiuVr = 



1 * 2 a-^) \ mm ). (8) 



2a 2 

\ 2a 2 2 ( iT To 1 1 



2 1 2 

Obviously, to a very good approximation, this matrix is of the "democratic" form in the 2,3 sector. 



IV. THE SEE-SAW NEUTRINO MASS MATRIX 

The see-saw neutrino mass matrix m v can now be obtained from Eq. (^) using (§) and (§). But before, we rescale 
the parameters a and Mo according to a — > -a, Mq — > yM^. One then gets the same expression for the light neutrino 
mass matrix m v as in the special case in which x and y in (Q) and (|^) are equal to 1: 



2 / —a 2 a 2 aa aa 

™»(Mo) = -t^tt a<J - 1 1 I (mt(M )f (9) 



2oM 



aa 1—1 



The basis chosen is still the basis in which mjj and m® zrac are diagonal matrices. Thus, we are not yet in a basis in 
which the charged lepton mass matrix is diagonal. We also have to use the renormalization group equation to go from 
the scale Mo down to the scale of the Z boson. But one can calculate from (^) the eigenvalues of the light neutrinos 
at Ma. As we will see below, the pattern of these eigenvalues will not change on the way down. 
The eigenvalues are 

(m t (M )) 2 fa 3 aa 4 

mi(M Q ) = - 
m 2 (M ) 



M \y/2 4 
(m t (M )) 2 fa 3 ac 



M V\/2 4 

ms( M 0) . (10) 

Mo a 

One can now identify the mass differences m|— mf rs m^—m 2 , with Am 2 observed in atmospheric neutrino experiments. 
The Super-Kamiokande result Am 2 sa 3 • 10 _3 (eV) 2 gives for aM^: 

aM « 5 • 10 11 GeV (11) 

For the ratio of the difference of squared masses, one gets 

5=4^1 = ^+0^) (12) 



2 2 

mi — m 



V2 
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Thus, the mass difference to 2 , — m\ relevant for solar neutrino oscillations turns out to be 

m\ - m\ w 4 • l(rV (eV) 2 . 



(13) 



Interestingly, we are left with the single parameter a only. 

As long as the off-diagonal elements of the charged lepton mass matrix are not much larger in magnitude than 
the off-diagonal elements of the down-quark mass matrix (see the next section), the diagonalization of (||) provides 
already an estimate for the neutrino mixing matrix U. The orthogonal matrix Om which diagonalizes m v {Mo) is 



M = 




(14) 



Because a/8 is a very small number, one has the remarkable result that the neutrino mixing matrix is of the bimaximal 
type |Q and very little_dependent on the mass difference of the 2 lightest neutrinos in the range acceptable for solar 
neutrino oscillations 



V. THE CHARGED LEPTON AND THE DOWN-QUARK MASS MATRICES 

We know the eigenvalues of the charged lepton mass matrix from experiment but not its form in a basis in which 
the Dirac neutrino and the up-quark mass matrices are diagonal. The down-quark mass matrix on the other hand is 
known to a good extent from the measured down-quark masses and the Cabibbo-Kobayashi-Maskawa mixing angles. 



In grand unified theories one obtains relations between to# and mp |10|, p5| , but the details depend on the Higgs 
representations and Yukawa couplings. The well-known relation obtained in SU(5) j|, |ll| 

m T E w m D (15) 

cannot be exact because of the different mass ratios in the quark and lepton sector. Below the scale Mo both wie and 
tod can be taken to be hermitian matrices. In the model considered here this should still be possible when reaching the 
scale Mo, which is much smaller than the unification scale where right-handed vector bosons are effective. Otherwise 
we would have to introduce more structure and thereby more unknown parameters. From the relations in grand 
unified theories such as jig ) one can then conclude that the mixing for charged leptons are similar to the mixing of 
quarks, i.e. not very largefx The mixing matrix U for neutrinos (at Mo) will then not much differ from Om given in 



To be more specific I will consider the hypothesis that at Mq, to^ and mu can be diagonalized simultaneously. In 
other words, in our model both matrices, if put into hermitian forms, commute as do the matrices m® lrac and mjj 



[m E , m D ] = . (16) 

On the one hand, this condition is weaker than (|l^) since the matrix elements of to|J and mo do not need to be equal 
or nearly equal. On the other hand, it is more strict since ( p^ ) allows a complete calculation of ras from the lepton 
masses and tod (or the CKM matrix). Of course all this has to be done at Mo and should then be scaled down by 
the renormalization group equations jl9| to the weak scale, i.e. the scale of the Z-boson. For mrj(Mo) we choose the 
matrix 

/ 0.7cr 3 il.54cr 2 ~ia 2 \ 
m D (M ) = -il.Ma 2 -a/3 i0.8a m b (M ) , (17) 
\ ia 2 -i0.8a 1 / 

which together with mu of Eq. (||) leads to mass eigenvalues at the weak scale within the experimental uncertainties 
(f|). It also gives at mz 



3 In this respect the present approach differs decisively from models which use asymmetric matrices by proposing for the quark 
sector large "right-handed" mixings of physical relevance correlated with a large mixing in the "left-handed" lepton sector wA , 
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|Kd|=0.22, \V cb \= 0.039, \V ub \= 0.0032 . (18) 

The phases chosen in ([Tt]) are such as to obtain "maximal" CP- violation Jl5|,|9| in the quark sector which may or 
may not be supported by ongoing experiments. 

An important matrix depending on the off-diagonal elements of mo and which is relevant for the amount of CP 
violation in the quark sector, is the commutator |20|] 

[mu, m D ) = ^C q {M ) . (19) 

The analog commutator for leptons is 

[m^ rac , m|] = -Ce(M ) . (20) 
i 

As an alternative to ( |l6| ) one can require that the off-diagonal elements of m E and have the same origin, i.e. arise 
from the same Higgs field: 

C e (M ) = ^^-C q (M ) . (21) 
m b (M ) 

The condition ( pl| ) , with x — 1 , mp from (|l^) , and the known charged lepton masses allows another calculation of 
m E {M ) 

-0.66cr 3 -il.Ma 2 ia 2 
m E (M ) = | il.54a 2 -a -i0.8a | m T (M Q ) (22) 
-ia 2 i0.8cr 1 

This matrix differs somewhat from the one obtained from ( |l6| ) but leads again to small charged lepton mixing angles. 

It would be interesting if another basis- independent relation between the mass matrices would exist, one which 
involves the lepton mass matrices only. Since Mr has no diagonal elements in our model one could speculatively 
assume 

^] = l Wt(M ° A ) ^ (M ° ) .M fl (Mo) • (23) 
i Mo 

Also this relation could be used to calculate the- The result, although somewhat different from the one obtained for 
m,E in (J22) fixes the phases of and leads again to charged lepton mixing angles of similar small magnitudes. 



VI. THE CALCULATION OF THE NEUTRINO MASSES AND THE NEUTRINO MIXING ANGLES 

At the scale Mq the mass matrix for light neutrinos is taken from Eq. (JsJ) . After the diagonalization of the mass 
matrix for charged leptons at this scale 

m E (M ) = U E {M )m% a9ona \M ) U^(M ) , (24) 

m u (M ) has to be transformed accordingly. The new neutrino mass matrix is then 

m„(M ) = (U E {M )) T m„(M ) U E {M ) . (25) 

This new neutrino matrix has now complex elements and will lead therefore to CP-violation effects in neutrino 
processes. m„(Mo) can now be scaled down to the weak scale by the corresponding renormalization group equation. 
The evolution equation according to the standard model is pl| : 

(Anf-^-mv = (-3.92 + 2A)m„ 

+ — 2 Tr[3mi/mjj + Snipm^ + ra E m E \m v 

-—^{fh v m E m E -\- (m E m E ) T m v ) . (26) 
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X = X(t) denotes the Higgs coupling constant related to the Higgs mass according to m 2 H — Xv 2 with v — 246 GeV. We 
take mii{mz) = 150 GeV for the numerical calculation. Solving (|26| ) allows to obtain the neutrino mass matrix m v 
at the scale of the standard model. The neutrino mixing matrix U = U(mz) can then be obtained by diagonalizing 
the hermitian matrix 77i*m y : 

m*(m z ) m u (m z ) = U(m* v m v ) dia9 ° nal tf . (27) 

The diagonal matrix 

^diagonal = jjT Vn„(m;f)Z7 (28) 

then provides the (complex) neutrino mass eigenvalues. By introducing the diagonal phase matrix 6 which consists 
of the phase factors of m diagonai divided by 2, U can be redefined U — > U<f>. The phase factors in (|2q ) cancel and the 
so obtained unitary matrix U expresses the neutrino states v e , v ^ , v T by the neutrino mass eigenstates according to 

Eq - ©■ 

The change of the mass matrices m\j , mo, tue and m v between mz and Mq depends, of course, on the renormal- 
ization group equations. To find the expression ( p~7| ) for mo (Mo) from the - approximate - knowledge of mo(rnz) I 
used the renormalization group matrix equations of the standard model [^9|. It leads to a change of \V c b\ by about 
10%. When using the renormalization group equation of the minimal supersymmetry model fl9[ |, the main difference 
to the standard model case lies in the overall factors mfc(M), m t (M), which can easily be adjusted. For m l/ (M) the 
formula (26) has to be replaced by the corresponding supersymmetry formula given by Babu et al. pl|. 



VII. RESULTS AND DISCUSSION 

It is straightforward to calculate from (p|),(^6|),(p^|),(p4])-(p8|) the light neutrino masses and the unitary mixing 
matrix U in terms of the parameters Mq and a. In a very direct way these two parameters determine the 2 mass 
differences in the 3 neutrino scenario. For the evolution of the mass matrices mjj ,mo,mE and m v as a function 
of the scale we applied the renormalization group equations according to the standard model and according to the 
minimal supersymmetry model. It turned out that for fixed Mq and a there is no noticeable change in the ratio of 
the 3 eigenvalues of m v when going from Mq down to mz- Even when the squared mass difference of the two lightest 
neutrinos is taken to be very small their mass ratio remains unchanged. Since both masses have opposite signs as 
seen in (|l0|), radiative corrections are ineffective as proved in ref. p2|] . 

Interestingly, it also turned out that the neutrino mixing matrix U is practically independent of the scale parameter. 
Moreover, the model predicts U to depend only very little on the difference of the neutrino masses. The reason is that 
in the mixing matrix (|lj) the parameter a plays a minor role only. Because also the mixing angles of the charged 
leptons are small as described in section 5, the matrix U of our model is of the bimaximal form |l4|] for all mass 
differences of interest in neutrino oscillation experiments^]. 

The deviations from bimaximal mixing are therefore almost exclusively due to the mixin gs o f charged leptons in 
the basis in which m® trac is diagonal. Using for mE{Mo) the result obtained from ( [l6|) and (|17|), one gets 

/ 0.69 0.71 0.15 \ 

Abs[U{m z )} = 0.51 0.51 0.69 (29) 
\ 0.51 0.49 0.71 / 

i.e. nearly pure bimaximal mixing. The amount of CP violation depends evidently on the corresponding CP-violating 
phase in the quark sector. With the phase choice taken in (|l7|), which at the weak scale provides an acceptable form 
of the conventional unitarity triangle for quarks, one obtains a neutrino unitarity triangle with angles: 

a„«77°, (3 U &17°, 7,-86° . (30) 

The form of this triangle is sensitive to the amount of CP violation in the quark sector. When using instead of ([l6]) the 
condition ( |2l| ) i.e. mE(Mo) from (^), (or from the condition (p3|)) the triangle becomes more flat (a v « 83°, (3 U ~ 
6°, j v « 91°) and \U e , 3 \ « 0.06. 



4 One may note here that for a small value of a, a « a, which describes a mass squared difference between the lightest neutrinos 
of order 10" 10 (eV) 2 , m„ in (^) takes a form not much different from the neutrino mass matrix suggested in Q and gives similar 
results. 
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In neutrino-less double (5 decay experiments p4 
predicts the small value 



one measures the neutrino mass averadge | < rh > e ^ e |. The model 



< m > f 



0.002 eV 



(31) 



As we have seen, the model presented here is predictive . Essentially, only the two mass differences between the 3 
neutrinos can be choosen freely. It predicts the bimaximal mixing form for the neutrino mixing matrix independently 
of the mass difference responsible for the solar neutrino oscillations. It is a minimal model because the parameter 
which describes the hierarchy in the quark sector could be used also for the lepton mass matrices, in particular, for 
the heavy neutrino mass matrix, providing in this way for an intimate connection between quarks and leptons. It is 
a minimal model also because only the simplest assumptions are made for the connection between m® lrac and mjj 
and between the and mp. The idea is that specific particles and interactions acting at the grand unification scale, 
besides producing the heavy neutrinos at Mo, may play a minor role at this lower scale. To give Mr the form of a Zee 
matrix is, of course, a strong assumption but may seem not unreasonable for dynamically generated singlet neutrino 
masses 4 to 5 orders of magnitudes below unification. One obtains for the heavy neutrinos two almost degenerate 
mass eigenstates which are strongly mixed. Accordingly, in the special basis in which Mr is diagonal, the matrices 
mu,mD,m® lrac and wle are of a "democratic" form in the 2,3 sector (see Eq. (||)). In other words, in this specific 
frame the Higgs field acts with (almost) equal strength on the members of the second and third family. However, it 
is worth pointing out that one can relax the Zee ansatz by including a 33 element to Mr. Because the inverse of Mr 
enters the see-saw formula even a very large 33 element (of order 10 2 ) still leads to near bimaximal mixing as long 
as the ratio (Mfl) 12 /(Mfl) 13 used in (g) is kept fixed. Thus, the oscillation properties of the light neutrinos cannot 
distinguish this strongly hierarchical mass matrix with small mixings of the heavy neutrinos^ from the mass matrix 
ofEq. @). 
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